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Eigenvalue cut-off in the cubic-quintic nonlinear Schrodinger equation 
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Using theoretical arguments, we prove the numerically well-known fact that the eigenvalues of all 
localized stationary solutions of the cubic-quintic 2D+1 nonlinear Schrodinger equation exhibit an 
upper cut-off value. The existence of the cut-off is inferred using Gagliardo-Nirenberg and Holder 
inequalities together with Pohozaev identities. We also show that, in the limit of eigenvalues close to 
zero, the eigenstates of the cubic-quintic nonlinear Schrodinger equation behave similarly to those 
of the cubic nonlinear Schrodinger equation. 
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I. INTRODUCTION 

The nonlinear Schrodinger equation (NLSE) has been 
widely used in modeling nonlinear wave dynamics in 
many physical scenarios, such as nonlinear optics ll|, 
2j, plasma physics Q, Bose-Einstein condensates |4j, 
biomolecular dynamics @ , and others [|| 0] • The sim- 
plest scalar NLSE is of the form 



dip 



= -M> + / 



i/). 



(i) 



where ip is a complex field denned usually on the whole 
R" with n = 1,2,3 and / describes the nonlinear re- 
sponse of the medium. Many different types of nonlinear- 
ities f(\ip\ 2 ) arise in the different physical fields of appli- 
cability of the equation including: power-law, saturable 
and nonlocal nonlinearities to cite a few examples. The 
most relevant one, for which many theoretical and analyt- 
ical studies of NLS equations have been done is the clas- 
sical cubic nonlinearity f{\ip\ 2 ) — g\ip\ 2 , both because of 
its direct interest and also because it corresponds to the 
simplest nonlinear response proportional to the square of 
the involved field (namely the light intensity in optics, 
the number of particles in BEC applications, etc.). 

One of the simplest extensions of the cubic NLS is the 
so-called cubic-quintic (CQNLS) model, which in nor- 
malized units, is 



i^ = -Ayj-g\yj\ 2 iP + h\ib\ 4 yj. 



(2) 



The CQNLS equation is another universal mathemati- 
cal model describing many situations of physical inter- 
est and approximating other more complicated ones. As 
examples it arises in plasma physics J9L |l0| , condensed 
matter physics |ll|. nuclear physics [Q, Bose-Einstein 
condensation [13j ]. etc., but probably the application of 
the model which has attracted more attention in the last 
years is description of the propagation of paraxial beams 



in certain nonlinear optical media. Many different op- 
tical materials have a refractive index that can be well 
described by a cubic-quintic nonlinearity such as some 
semiconductors and doped glasses (e.g. AlGaAs [14j and 
CdS^Sei-a, (lal. the polydiacetylene para-toluene sul- 
fonate (PTS) [16J, chalcogenide glasses [13], some trans- 
parent organic materials [18| or even media with complex 
susceptibilities induced by Electromagnetically Induced 
Transparency [191 ]. 

We will consider localized stationary solutions (i.e. 
solitary waves or solitons) of Eq. ©. Thus, taking 
i()(x, t) — u(x)e ll3t we will study solutions of 



(3u = Au + (g\u\ 2 - h\u\ A ) ■ 



(3) 
and vanishing 






with u a complex function defined on R 2 
at infinity. 

The stationary solutions of Eq. ([3]) and their stabilit 
properties, have been studied in many papers 
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A well-known fact is that no solitary waves are found 
beyond a maximum value of the eigenvalue /3 = j3*. In 
fact for (3 < /3* one finds many solutions with different 
widths corresponding to very different values of the norm 
of the solution Np = J R2 \up\ 2 dx. In fact, when j3 — ► /3*, 
N/3 — ► oo. Thus, there is a cut-off in the permitted eigen- 
value of localized solutions of the CQNLS. The existence 
of this cutoff and its specific value as a function of the pa- 
rameters was discussed in many papers. For fundamental 
states it was studied numerically [37] and by approximate 
variational methods with a supergaussian ansatz in Ref. 
[231 ] and also in Ref. [28| • For vortex states the same 
problem was studied by means of variational [331 ] . nu- 
merical [33[ and analytical methods [26|, [3lJ • 

In this paper we provide a theoretical support for 
previous works dealing with the problem of the cutoff 
[23J, [26|, [28|, [3l|, [33j and prove in a rigorous yet simple 
way that localized stationary states of the cubic-quintic 
NLS equation exist only for eigenvalues on a finite inter- 



val. The obtained result holds for all type of localized 
solutions, e.g. ground, vortex or dipole states, of the 
cubic-quintic NLS equation. 

The article is organized as follows: first in Section |TT] 
we introduce the NLS equation to be considered, the in- 
tegral quantities and the inequalities to be used in the 
subsequent sections. Section IIIII contains the derivation 
of the upper limit for the eigenvalues and analysis of the 
behavior near cut-off. Next, in Section ITVl we describe 
the related problem of the asymptotic behavior of the 
stationary states in the NLSE limit. Finally, in Sec. [V] 
we summarize our results. 



II. STATEMENT OF THE PROBLEM 

The typical situation for the cubic-quintic nonlinear- 
ity corresponds to the case where we have a combination 
of focusing cubic and defocusing quintic term. The in- 
terplay of focusing and defocu sing nonlinearity in that 
case prevents the wave collapse [39] and is responsible for 
the liquid-like features of the localized stationary states 

Hi El pi . 

In this section we will consider the case of two spatial 
dimensions which is the one relevant to nonlinear optics 
and many other applications of the model. 

The equation we are dealing with can be written as 



This relation can be rewritten in terms of X = N/K 



(3u — Am + gu 3 - hu 5 , 
with [3, g,h > 0. Our analysis will use the energy 

K = -0N + gU 2 - hU 3 
and Pohozaev identities [JO . |41| 

= pN-£u a + ~U a , 

where the moments of u are defined by: 



K 



cTrlVd 2 , 



/' 



N = ld n ru 2 , 



Ui --■■ ld n ru 2 \ j = 2,3. 



(4) 
(5) 

(6) 

(7a) 
(7b) 
(7c) 



III. CUTOFF IN THE CUBIC-QUINTIC NLSE 

A. Existence of cut-off 

By direct calculations from |5j) and ([6]) we get the fol- 
lowing bilinear relation between the moments 

AhK 2 + (3. 9 2 - 16/%) NK + [3 (l6/3/i - 3g 2 ) N 2 

+ g 2 h(NU 3 -U 2 )=0. (8) 



/3(/3-/3*)A 2 + (/3» 



/3)x + ^(i + r) = o, 



where 



g 2 NU 3 - U 2 
4 K 2 


/3* = 


_ 3 9 2 
16 h 



(9) 



(10) 



It should be noted that it follows from the Holder in- 
equality 

NU 3 > U 2 , (11) 

that r > 0. Now, since our quadratic equation ([9]) must 
have real roots, after calculating its discriminant 



P = (/3* 
one can conclude that 



P) W* + /ST) 



/3</3* 



(12) 



(13) 



Thus, there is an upper bound on the eigenvalues, 
3g 2 /16h and localized stationary states of Eq. ((4]) ex- 
ist only when (3 lies within the interval 



-'^ 



(14) 



This result agrees with previous numerical and approx- 
imate calculations for this quantity [23|, [26|, [28|, [Ml [331 
but here is obtained rigorously using simple arguments. 
In fact, the comparison with the numerical results of Rcf. 
[231 . |28L 1331 ] shows that our bound is optimal and matches 
closely the numerics. 



B. Behavior near cut-off 

A natural question that arises, after establishing the 
bound p3[) is what occurs with solution when the pa- 
rameter (3 approaches the cut-off limit. The aim of this 
section is to obtain a bound for the quantity N which 
allows us to show that iV — > +oo when j3 — ► /3*. 

Solving Eq. © we obtain that 



N _ 1 
K ~ 2/3 




/3«+I7? 



(15) 



The sign before the radical is determined by the identity 



2[3N-K= X -gU 2 > 0, 



(16) 



which is a direct consequence of Eqs. ([5]) and ©. Since 
r > 0, Eq. (fl5|) leads to the inequality 



/3* 



N 1 

K > 2/3 V #*-/3' 



(17) 



On the other hand we can obtain another type of bounds 
for K and A. From Eqs. ([5]) and ^ we get 



where 



~U 3 = /3N-K, 

K = ~gU 2 - hu 3 . 

Combining these relations and the identity 
dK 



dft 



N, 



(18a) 

(18b) 

(19) 



which can be derived by differentiating Eq. (fj[]) and ap- 
plying again Eqs. © and ©, we conclude that K is a 
monotone function. These leads to the fact that 



K = when ft = 







(20) 



and respectively K = J N dft. Now let us recall the par- 



ticular form of Gagliardo-Nirenberg |43j, |46j, |47j inequal- 
ity which in our case can be written as 



U 2 < C GN KN, 



(21) 



where C GN is the optimal constant for the Gagliardo- 
Nirenberg inequality in two dimensions. Applying this 
inequality to Eq. (|18b[) we get 



A> 



Combining Eqs. (|22|) and (fT9l) we get the inequality 

2 



(22) 



K> 



Cgn<? 



ft. 



(23) 



Substituting this estimate into (fT?| we obtain finally 



N > 



ft* 



C G 



ft* -ft' 



(24) 



which guarantees that A — ► +oo as ft — > /3* = 
3<? 2 /(16/i). Thus we can see that the cut-off phenomenon 
is a manifestation of the blow-up of the norm (note that 
A is nothing but the L 2 -norm of u). 

Another consequence of the estimate (|24|) is the fact 
that A is bounded from below by \jC GN when ft — » 0. 



IV. THE LIMIT /3 -> 

It is posible to analyze the asymptotic behavior of so- 
lutions of Eq. ((4| when ft — ► 0. To do it we use the 
rescaling symmetry of our equation, specifically we will 
use the fact that if u is a solution of Eq. (fj|]) then the 
function u given by u(r) = £,u(r]r) solves 



, i , ° 



/3u = Am + .<?u J - /iu 



(25) 



V 



'I 



ft = V ft, 9= T2-9, h= -j/i 



(26) 



Calculating the moments of the function u we arrive at 
the following rescaling relation for K and N treated as 
functions of the parameters ft, g and h: 

K(ft,g,h) - 1^ (y^ |_<? ; |_ A (27a) 

2 / 2 2 \ 

N(ft,g,h) = ^ A (y/3, ^, l^/ij , (27b) 
Further, choosing 



2 1 t i 9 

v =-ft> e = -ft 



we can obtain 



9 \ 9 



(28a) 



A(/3, ff ,/i) - ijv(l,l,^V (28b) 

Considering the limit ft — » we find the following asymp- 
totic formulae for N[u] and K[u]: 



N[u] = -N[<P}+0{ft) 
9 



(29a) 



K[u] = ^K[0] + ^[/ 3 [0] + O(/? 3 ), (29b) 

9 Sg* 

where <j) is the Townes soliton, i.e. the solution of 

A0 = 0-0 3 . (30) 

These results make it possible to deduce that in the 
asymptotic region as ft — > the behavior of solu- 
tions of cubic-quintic Schrodinger equation is similar to 
the behavior of eigensolutions of the cubic Nonlinear 
Schrodinger equation. 



V. CONCLUSIONS 

Using general mathematical arguments we have shown 
that the localized stationary solutions of the cubic- 
quintic Nonlinear Schrodinger equation exhibit an eigen- 
value cut-off. This result holds for all type of localized 
stationary solutions, e.g. fundamental, vortex or dipole 
states. The cut-off is in excellent agreement with previ- 
ous approximate or numerical results. We have also ob- 
tained a lower bound for the number of particles in the 
eigenstates that diverges exactly at the cutoff ft = ft*. 
Finally, in the limit of small eigenvalues we obtain the 
result that the eigenstates of the cubic quintic nonlinear 
Schrodinger equation behave in the same way as those of 
the cubic nonlinear Schrodinger equation. 

Our results complement present knowledge on one of 
the key models of mathematical physics and support in a 
rigorous yet simple way previous numerical observations. 
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